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We consider the leading twist T-odd contributions as the dominant source of the cos 2<f> azimuthal asymmetry in 
unpolarized pp — > fi/i + X dilepton production in Drell-Yan Scattering. This asymmetry contains information on 
the distribution of quark transverse spin in an unpolarized proton. In a parton-spectator framework we estimate 
these asymmetries at 50 GeV center of mass energy. This azimuthal asymmetry is interesting in light of proposed 
experiments at GSI, where an anti-proton beam is ideal for studying the transversity properties of quarks due to 
the dominance of valence quark effects. 



1. Introduction 

One of the persistent challenges confronting 
the QCD parton model is to provide a theoret- 
ical basis to understand the experimentally sig- 
nificant azimuthal and transverse spin asymme- 
tries that emerge in exclusive, inclusive and semi- 
inclusive processes [1I2I3I4I5I5] . Generally speak- 
ing, the spin dependent amplitudes for the scat- 
tering will contribute to non-zero transverse sin- 
gle spin asymmetries (SSA) if there are imaginary 
parts of bilinear products of those amplitudes 
that have overall helicity change. For two-body 
exclusive reactions, SSA requires there to be an 
imaginary part of the product of an helicity non- 
flip with an helicity flip amplitude. For inclusive 
reactions, the same conclusion can be reached by 
taking the amplitudes as two-body helicity ampli- 
tudes for the production of a fixed hadron and a 
state \X >. Through the generalized optical the- 
orem, SSA in inclusive reactions can be related to 
discontinuities in helicity flip three-body forward 
scattering amplitudes |7j . That is essentially spin 
kinematics. Dynamically there must be quantum 
field theory contributions to the relevant ampli- 
tudes. In perturbative QCD (PQCD), applicable 
to the hard scattering region, to obtain an imagi- 
nary contribution to quark and/or gluon scatter- 
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ing processes demands introducing higher order 
corrections to tree level processes. One approach 
incorporates the requisite phases through inter- 
ference of tree level and one-loop contributions 
in PQCD in an attempt to explain spin asymme- 
try in A production [Sj. On general grounds the 
helicity conservation property of massless QCD 
predicts that such contributions are small, going 
like a s m/Q, where a s is the strong coupling, m 
represents a non-zero quark mass and Q repre- 
sents the hard QCD scale |9I8| . Such contribu- 
tions have failed to account for the large SSA ob- 
served in A production 1 . On the other hand, 
the twist three quark-quark and quark-gluon cor- 
relations described in [TOj and [ITj hold promise 
to describe the phenomena at large px- 

However, considering the soft contributions to 
hadronic processes opens up the possibility that 
there are non-trivial transversity parton distri- 
butions that can contribute to transverse spin 
asymmetries. Ralston and Soper introduced J^l 
the now well known chiral-odd transversity trans- 
fer distribution function |lrill4| h\{x) which can 
play a role in doubly polarized Drell-Yan pro- 
cesses |15) . h\{x) can also be measured in semi- 
inclusive deep inelastic scattering (SIDIS) [16) . 
For SSA in SIDIS transverse momentum must be 
acquired to lead to appropriate helicity changes. 
In describing transverse asymmetries this is par- 
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ticularly important when the transverse momen- 
tum is sensitive to intrinsic quark momenta. Here 
the effects are associated with non-perturbative 
transverse momentum distribution functions |17) 
(TMD), where transverse single spin asymme- 
tries indicate so called T-odd correlations be- 
tween transverse spin and longitudinal and in- 
trinsic quark transverse momentum. The T-odd 
distributions |18I19I2L)| are of importance as they 
possess both transversity properties and the nec- 
essary phases to account for SSA and azimuthal 
asymmetries |16I21I22) . Formally, these phases 
can be generated from the gauge invariant def- 
initions of the T-odd quark distribution func- 
tions |23I24I25I2^| . Asymmetries that involve T- 
odd TMD and fragmentation functions are in- 
dicative of a rich set of correlations among trans- 
verse momenta of quarks and/or hadrons and the 
transverse spin of the reacting hadrons and/or 
quarks. In contrast to the SSA generated from 
the interference of tree-level and one loop correc- 
tion in PQCD, such effects go like a s <k±>/M, 
where now M plays the role of the chiral sym- 
metry breaking scale and k± is characteristic of 
quark intrinsic motion. 

T-odd distributions only exist by virtue of non- 
zero parton transverse momenta |19I27I20] . They 
also correspond to distributions that would van- 
ish at tree level in any T-conserving model of 
hadrons and quarks. In this sense they are simi- 
lar to the decay amplitudes for hadrons that in- 
volve single spin asymmetries which are non-zero 
due to final (and/or initial) state strong interac- 
tions. Their existence was suggested by Sivers 
to account for the significant SSA in inclusive re- 
actions (e.g. pp1 — > irX) |18I19) . by Collins in 
SIDIS HSI, and by Boer [2U in Drell-Yan scatter- 
ing. 

A great deal of progress has been made in 
the last several years, following the realization 
that this link goes further, as Brodsky, Hwang 
and Schmidt first showed [22| and several re- 
searchers generalized thereafter |24I25| . Detailed 
model calculations of these functions have been 
performed in parton inspired spectator-models of 
quark-hadron interactions where absorptive ef- 
fects are generated by interference between gluon 
loop corrections to tree-level TMDs. These loop 



calculations were applied to SIDIS |22l24l28l25l 
I29I3()I31I32| and Drell-Yan processes |33I34I35| . 

thereby giving rise to predictions for SSA and az- 
imuthal asymmetries. 

The importance of quark distribution func- 
tions was recognized some 35 years ago by Drell 
and Yan |3fi) when they considered high energy 
hadron scattering that produces large invariant 
mass lepton pairs as a fundamental probe of 
quark-antiquark distribution functions. Further- 
more, considering various asymmetries and po- 
larization phenomena in Drell-Yan processes can 
uncover relevant products of spin dependent dis- 
tributions j2J. Indeed, Drell Yan pp scattering 
is a preferred reaction to study the the role that 
T-odd quark distribution functions play in the 
transverse spin structure of the proton through 
spin and azimuthal asymmetries in QCD |21j . 
That is the direction we pursue herein. 

2. Drell-Yan and T-Odd Correlations 

At the parton level the Drell-Yan cross section 
will receive contributions from quark-antiquark 
annihilation into the heavy photon. In unpo- 
larized Drell-Yan scattering early cross section 
data as a function of the transverse momentum 
of the muon pair indicated deviations from the 
Bjorken scaling prediction |37I38| . The impli- 
cation was that the collinear approximation was 
insufficient to describe the data |39I40| . Trans- 
verse momentum of a parton arises due to hard 
Brcmsstrahlung of gluons, which is calculable 
from PQCD when the momentum transfers are 
large 41 . On the other hand, quark confinement 
implies that quarks have soft, primordial or in- 
trinsic transverse momenta k±. This latter ef- 
fect is significant at low transverse momentum, 
q T <^ Q. q T dependence has been incorporated 
into the factorized Drell-Yan model |12I42 | by ex- 
tending the parton probability distribution to be 
a function of k± ^7] 3 

J dk x V(k x ,x) = f(x). (1) 

If the parton distributions within the incoming 
hadrons have transverse momentum dependence, 

3 New work on the factorization theorem for Drell-Yan can 
be found in 1431 . 
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there will be a continuum of values of their k±_ for 
which a time-like photon of fixed 4-momentum 
will be formed. Ignoring or summing over spin 
(and the lepton pair orientation), the k± depen- 
dent distribution functions appear in the differ- 
ential cross section, 



da 



Ana- 



^e 2 



d 2 k. 



d 2 P A 



dq 2 dyd 2 q T " 3Q 4 ^ a , 

(k ± +p ± - q T ) f a/A (x,k±)fa/B{x,P±) 

where f a /A{xi, kj_) is a distribution function for 
a quark a to be found in hadron A with trans- 
verse momentum fcj_ and longitudinal momentum 
fraction x\ and / is the corresponding anti-quark 
distribution in hadron B. 

Once transverse momentum dependence of par- 
ton distributions enters the picture of scattering 
processes a much larger set of transverse mo- 
mentum distribution (TMD) and fragmentation 
functions 27] become possible and relevant, par- 
ticularly for spin asymmetries. Among such func- 
tions are the possible leading twist T-odd quark 
distribution |18| and fragmentation functions . 

In SIDIS with an unpolarized target, an expec- 
tation value of isx ■ (P x fe_i_), indicates a T-odd 
correlation of transverse quark polarization with 
the proton's momentum and the intrinsic quark 
transverse momentum in an unpolarized nucleon 
while isx ■ (p x Ph±), corresponds to that of a 
fragmenting quark's polarization with quark and 
transverse pion momentum, PhX- These correla- 
tions enter the unpolarized cross-section through 
convolutions with J2U| and the Collins frag- 
mentation function H^- \^ . The resulting cos 2(f) 
asymmetry is not suppressed by 1/Q, where Q 
represents the scale in SIDIS - the T-odd contri- 
bution is at leading twist |2UI29I3L)| . An analo- 
gous unpolarized double T-odd azimuthal asym- 
metry enters the Drell-Yan process For the 
Drell-Yan process the angular dependence can be 
expressed as 

-l 



dN 

In 



da 



dQ 2 dydq\, 



da 



dQ 2 dydq\d£L 



1 



4tt A + 3 



1 + A cos 2 9 + /i sin 6 cos <j) 

+ ^sin 2 6>cos2</>). (2) 



The solid angle Q refers to the lepton pair orien- 
tation in the pair rest frame relative to the boost 
direction and the incoming hadrons' plane |42) . 
A, fx, v are functions that depend on s, x, tti 2 ^, q T , 
the center of mass energy, the fraction of quark 
momentum in the hadron , the invariant mass of 
the produced lepton pair, and the transverse mo- 
mentum of the dimuon pair. 4 

By keeping the Q dependence in the convolu- 
tion, Collins and Soper could project out the 
photon angular dependence in Eq. and obtain, 
among the other asymmetry functions, a (spin av- 
eraged) T-even contribution to the cos2</> asym- 
metry 



ZJI Ea [WA fl (X, fcjj/l (X,P±) 

Ea e l T k ±-)fl( x ,P±)) 



(3) 



where W4 = 2(h ■ (k± — p±)) 2 — (fej_ — P±) 2 and, 
h = q T /Qr- Indeed the earliest theoretical ex- 
planation for azimuthal asymmetries was given 
by Collins and Soper 's estimate of v in Eq. ©, 
which is a kinematic non-leading twist contribu- 
tion. However, all of the asymmetry functions, 
/i, A and v 1 were shown to have parton model 
contributions. Taking into account NLO and 
NNLO 03] the QCD improved parton model pre- 
dicts 1 — A — 2v = 0, the so called Lam- Tung 
relation |45| . However, experimental measure- 
ments of Tip — > n. + fj,~X discovered unexpectedly 
large values of these asymmetries [213] compared 
to parton-model expectations resulting in a se- 
rious violation of this relation. Attempts to ac- 
count for the violation in terms of higher twist 
effects [46147] have been unsuccessful. More re- 
cently, Boer |2I] proposed that there is a domi- 
nant leading twist contribution to v coming from 
the T-odd transversity distributions hj^(x, k±) for 
both hadrons which dominates in the kinematic 
range, qx <C Q- The cos 2(f> azimuthal asymmetry 
in unpolarized pp — > /i + pT X would involve the 
convolution of the leading twist T-odd function, 
hi |21I33I34I3S| 



V-2 



E a e\T [w2 hj(x, fcx)fej L (^,P±)/(MiM 2 )] 

Ea e l T [hi x > kx)fl{x,Px)] 



4 We are working in the Collins-Soper frame where q T re- 
tains its meaning. 
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where W2 = {2h ■ k± • h ■ p ± — p ± ■ k±) is the 
weight in the convolution integral, T |21| . A 
simple model for these distributions, inspired by 
Collins' ansatz for the transversity fragmentation 
function led to Qt dependent v which could be 
fit to the low values of the -np data. Further work 
along those lines 33 34 35 incorporated a more 
realistic model for the T-odd functions, as first de- 
veloped in SIDIS [22] for the functions f^(x, k±) 
which were related to the h^(x,k±) in ref. [25] , 
The results were presented for pp scattering. 5 
This azimuthal asymmetry is interesting in light 
of proposed experiments at Darmstadt GSI |5(J| . 
where an anti-proton beam is ideal for studying 
the transversity property of quarks due to the 
dominance of valence quark effects |49| . Herein 
we extend our calculations for T-odd contribu- 
tions to the unpolarized Drell-Yan Scattering first 
reported in We perform a detailed analysis 
displaying qx and for the first time x, xp, and 
q (or m w ) dependence of this effect. In addi- 
tion we compare the double T-odd contribution 
to the conventional subleading twist T-even con- 
tribution [1*2] . 

3. T-odd transversity distribution 

For our purposes we consider projected 
from the correlation function for the TMD func- 
tion $(fc,P), 

+ ...} (5) 

that is 



J dk-Tr ((7^+75$) 



2£+ -^ k±j hi(x,k ± ). 



In our work on SIDIS we used a parton model 
within the quark-diquark spectator framework to 
model the quark-hadron interactions that enter 
$(fc,P) |25I29| and contribute to T-odd terms 
in the projection. The basic diagram, indicating 
the one loop gluon exchange and the eikonalizcd 

5 Very recently instanton induced effects have been inves- 
tigated ESI . 




+ 



H.C. 



H.C.. 



Figure 1. Above: Feynman diagram representing fi- 
nal state interactions giving rise to T-odd contribu- 
tion to Drell-Yan Scattering. Below: Quark-target 
scattering amplitude depicting the T-odd contribu- 
tion to the quark distribution function in the eikonal 
approximation. 



struck quark line arising from the gauge link is 
indicated in Figure pj Noting that parton intrin- 
sic transverse momentum yielded a natural regu- 
larization for the moments of these distributions, 
we incorporated a Gaussian from factor into our 
model |29I3(J| . The result was 



hi(x,kx) =Na s M 



(1 -x)(m + xM) 



n h (k 



2 . 

J.! ' 



(6) 



where 1Z is the regularization function 



— o vl -,- 2h ( fe i- A (°)) 



TZ(k ± ; x) = exp 

x (r(0, 26A(0)) - r(0, 26A(*i))) . 

J\f is a normalization factor determined with re- 
spect to the unpolarized u-quark distribution, ob- 
tained from the zcroth moment of 



fi(x,kx) = 



M{l-x) [(jn + xMf + k 2 ^ 



normalized with respect to valence distributions 
in ISU- Tl s {k\-x) = exp- 2bfc i. 

The asymmetry in the pp Drell-Yan process in- 
volves the convolution of the product of two T- 
odd distributions. The contribution to the double 
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T-odd azimuthal cos 2(f> asymmetry, v in Eq. Q , 
in terms of initial and final (ISI/FSI) state in- 
teractions of active or "struck" , and fragmenting 
quark |29I3U| is depicted in Figure ^ We 




Figure 2. v plotted as a function of qr for s = 
50 GeV 2 , x in the range 0.2 — 1.0, and q ranging 
from 3- 6 GeV/c. 



have numerically performed the convolution inte- 
grals and obtained values of the asymmetry v as 
a function of the variables, x, xp, qp, and q (or 

Before evaluating the convolution, the Drell- 
Yan kinematics demand some special attention. 
With x and x being the fractional longitudinal 
momenta of the quark and antiquark, there are 
some constraints: 

xx = t = Q 2 /s , — - — = r) = xp/2 and 
x = n + \f rj 2 + t 2 , x — — f] + \/r\ 2 + r 2 . (7) 

Due to the constraint xx — q 2 /s the allowed 
range of x is restricted for each q value, from 
Xmin — q 2 / s to 1. Furthermore, evaluating the 
convolutions of hj-hj- and /i/i for a sampling of 
x will not treat the X and the corresponding an- 
tiparticle structure functions symmetrically. So it 
is more appropriate to use the symmetrical vari- 
able, Feynman-x, xp = x — i or x = \{xp + 
y/x 2 F + 4q 2 /s) and x = \{-x F + y/x 2 F + 4q 2 /s). 
However, the allowed range of xp depends on q 




0.1 - 



Figure 3. v plotted as a function of q = m M(J for 
s — 50 GeV 2 , x in the range 0.2 — 1.0, and qr ranging 
from 0-3 GeV/c. 



(from -l/2(l-q 2 /s) to +1/2(1 -g 2 /s))-the vari- 
ables xp and q are not orthogonal. Since we aim 
to present partially integrated values of ap- 
proximating experimentalists' measurements, it is 
advantageous to work with orthogonal variables. 
We choose the variable 

with range from — i to +h, independent of q. Val- 
ues of the asymmetry fill the rectangular space of 
variables, C, q, qp- We have been careful with this 
choice because our model predictions have con- 
siderable structure in all 3 variables. Hence the 
meaning of a graph of v{q) or v(x) has particu- 
lar significance when comparing to experimental 
data. 

A crucial point in selecting these variables 
involves how experimenters determine various 
asymmetries and angular dependences, in order 
to maximize statistics when extracting possibly 
small effects like v. Events appear distributed 
over allowed regions (modified by experimental 
acceptances) of all three variables along with the 
[i pair angular variables, of course. To obtain the 
dependence on one variable, large bins are de- 
fined and event numbers averaged over those bins. 
How arc those results to be compared with the- 
oretical predictions |52| ? The two experiments 
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Figure 4. Contours of constant a; as a function of £ 
and q which ranges from to 3 GeV/c. 



for which data have been published have differ- 
ent ranges of variables \2l.'>\ . The binning proce- 
dures are not easily compared. To be most gen- 
eral and adaptable for future experimental com- 
parisons we have determined the value of v as a 
function of (,q,qr- We then integrate over pairs 
of those variables for particular ranges of the vari- 
ables. At s — 50 GeV 2 we take qr < 3 GeV/c and 
3 GeV/c < q < 6GeV/c, while £ always varies 
from -1/2 to +1/2. 

For v{qT) and v(q) the resulting values are 
shown in Figures [21 and [3J There is a correspond- 
ing shown in Figure [SJ To connect with the 
x or xf dependence we have to take the q depen- 
dence into account. For any qr the fixed xp val- 
ues form contours in the C, <7 plane (see Figure^}. 
So an integral of v over q for a fixed xf follows 
the relevant contour in Q, q and has a limit on 
t = q 2 /s of (1 — T^p). This limits the range of the 
q integral until the limiting value of the range at 
q = 6 GeV/c for s = 50 GeV 2 . Similarly, for any 
qr the fixed x values form asymmetrical contours 
in the £ ; 9 plane as shown in Figured The limit 
on t for a fixed x will be x( ^^ 1 ). The resulting 
values of v{x) are shown in Figure [HI Finally, in 
Figure[7|we plot the leading twist contribution to 
v as a function of Xf- 




Figure 5. u plotted as a function of C, for s = 
50 GeV 2 , qr ranging from 3 to 6 GeV/c and q from 
to 3 GeV/c. 



4. T-Even Contribution 

Long before the realization that there is a lead- 
ing twist 2 contribution to the Drell-Yan az- 
imuthal asymmetry, it was proposed by Collins 
and Soper |42| that the spin independent, trans- 
verse momentum dependent distributions /i and 
fi could contribute via Eq. It is important 
to compare this kinematic twist 4 contribution 
to the leading twist contribution (Eq. (0J shown 
above. We combined both convolutions to de- 
termine the magnitude of the shift. The addi- 
tional contribution for s = 50 GeV 2 to each of 
the partially integrated functions v is shown in 
Figures I2I3I5IKI as slightly higher curves. At most 
the additional contribution is around 3-4%. For 
higher s values the effect is even smaller, as ex- 
pected El]. 

5. Conclusion 

A perusal of the figures shows that the cos 20 
azimuthal asymmetry v is not small at center of 
mass energies of 50 GeV 2 . We estimated the lead- 
ing twist 2 and twist 4 contributions [32] . In Fig- 
ure [21 the T-odd portion contributes about 30% 
with an additional 3% from the sub-leading T- 
even piece. The distinction between the leading 
order T-odd and sub-leading order T-even con- 
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Figure 6. v plotted as a function of a; for s = 50 GeV 2 
qr ranging from 3 to 6 GeV/c and q from to 3 
GeV/c. 



tributions diminish at center of mass energy of 
s = 500 GeV 2 HH|. In Figure El v is plotted ver- 
sus x at s = 50 GeV 2 , where qx ranges from 2 
to 4 GeV. Again the higher twist contribution is 
small. 

Thus, aside from the competing T-even ef- 
fect, the experimental observation of a strong 
x-dependence would indicate the presence of T- 
odd structures in unpolarized Drell-Yan scatter- 
ing, implying that novel transversity properties 
of the nucleon can be accessed without invoking 
beam or target polarization. 

It should be noted that at order a s a complete 
analysis for the full range of qr would entail in- 
cluding gluon bremsstrahlung contributions |41) . 
Furthermore, collinear Sudakhov corrections have 
not been accounted for here [53] • A thorough 
explication of Drell-Yan dynamics would require 
more care with regions in which divergent contri- 
butions become important to address. For this 
study, however, we have considered the implica- 
tions of our model, unencumbered by subtleties 
at the edges of the phase space on which we con- 
centrate. 

We conclude that T-odd correlations of intrin- 
sic transverse quark momentum and transverse 
spin of quarks are intimately connected with stud- 
ies of the cos 2(f> azimuthal asymmetries in pp- 
Drell-Yan scattering. Due to the dominance of 




Figure 7. The leading twist contribution to v plotted 
as a function of xf for s — 50 GeV 2 qr ranging from 
3 to 6 GeV/c and q from to 3 GeV/c. 



valence quark effects we estimate that the pro- 
posed proton anti-proton experiments at GSI |50| 
provide an excellent opportunity to study the role 
that T-odd correlations play in characterizing in- 
trinsic transverse spin effects within the proton. 
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